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ABSTRACT 



We present the results of an 850 /im JCMT/SCUBA survey for dust around 
13 nearby solar mass stars. The dust mass sensitivity ranged from 5 x 10^^ to 
^ ■ 0.16 M®. Three sources were detected in the survey, one of which (HD 107146) 

(^Q ■ has been previously reported. One of the other two submillimeter sources, 

^ ■ HD 104860, was not detected by IRAS and is surrounded by a cold, massive dust 

O ■ disk with a dust temperature and mass of Tdust = 33 K and M^ust = 0.16 M^. 

■ The third source, HD 8907, was detected by IRAS and ISO at 60 — 87 /im, and 

has a dust temperature and mass of Tdust = 48 K and Mdust = 0.036 M®. We find 
that the deduced masses and radii of the dust disks in our sample are roughly 
^ . consistent with models for the coUisional evolution of planetesimal disks with 

embedded planets. We also searched for residual gas in two of the three systems 



with detected submillimeter excesses and place limits on the mass of gas residing 
in these systems. 

When the properties measured for the detected excess sources are combined 
with the larger population of submillimeter excess sources from the literature, we 
find strong evidence that the mass in small grains declines significantly on a ~ 200 
Myr timescale, approximately inversely with age. However, we also find that the 
characteristic dust radii of the population, obtained from the dust temperature 
of the excess and assuming blackbody grains, is uncorrelated with age. This is in 
contrast to self-stirred collisional models for debris disk evolution which predict 
a trend of radius increasing with age tage oc R^. The lack of agreement suggests 
that processes beyond self-stirring, such as giant planet formation, play a role in 
the evolutionary histories of planetesimal disks. 



Subject headings: circumstellar matter — planetary systems: protoplanetary 
disks — planetary systems: formation — submillimeter 
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1. Introduction 

The existence of dust around stars well beyond the protostellar phase is thought to arise 
from the collisions of planetesimals, the building blocks of planets, and is therefore of great 
interest for the clues they may reveal about the processes and timescales involved in planet 
formation (Backman & Paresce 1993). Indeed dusty disks are thought to be one of the most 
readily observable signatures of extrasolar planetary systems. Abundant dust is expected 
to be produced early on as part of the planet formation process itself, as planetesimals 
collide and grow into planetary mass objects (e.g., Kenyon & Bromley 2004a, b). Further 
dust production is expected as residual planetesimal disks grind themselves down to the 
low masses inferred for the Kuiper Belt today (e.g., Farinella, David, & Stern 2000). Giant 
planets in the system may sculpt the resulting debris, providing an indirect signature of their 
presence (e.g., Liou & Zook 1999; Ozernoy et al. 2000; Quillen & Thorndike 2002; Moro- 
Martin & Malhotra 2002). In the context of our solar system, the age range 10-100 Myr 
is of interest, since it overlaps the epoch of giant planet formation and includes the epoch 
of formation of large bodies in the outer solar system. Also during this time period, the 
interaction of giant planets with a residual planetesimal disk may lead to planetary orbital 
migration and sculpting of the planetesimal disk. 

Space based far-infrared telescopes have discovered and begun to characterize the prop- 
erties of dusty disks around main sequence stars (Aumann 1985, Spangler 2000) and it is 
one of the primary goals of the recently launched Spitzer Space Telescope. Ground-based 
sub-millimeter observations provide useful complementary information about the properties 
of such disks. First, the dust mass is readily measured because the emission is almost cer- 
tainly optically thin at long wavelengths. Second, by anchoring the Rayleigh- Jeans side of 
the spectral energy distribution (SED), contraints can be placed on the temperature and size 
of the emitting grains. Finally, the relatively high resolution that can be attained at these 
wavelengths allows investigations of the morphology of the emitting dust. Dust morphologies 
have been used to argue for the presence of perturbing planetary companions (Holland et al. 
1998, Greaves et al. 1998, Wilner et al. 2002). 

In this paper we present the results of a sub-millimeter search for dust debris surrounding 

13 solar mass stars that are targets of the Spitzer Legacy Science Project, "Formation and 
Evolution of Planetary Systems" (FEPS; Meyer et al. 2004). Since a goal of the FEPS 
project is to place our solar system in context, the parent sample covers spectral types 
corresponding to approximately solar mass stars spanning a range of ages from 3 Myr to 3 
Gyr. The submillimeter targets were selected from the FEPS sample with an emphasis on 
sources that are young (ages 10 to a few 100 Myr) and nearby (< 50 pc) based on estimated 
stellar ages for the FEPS sample (Hillenbrand et al. 2005) and Hipparcos distances. We 
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emphasized nearby sources in order to enable the detection of faint excesses or to place 
meaningful upper limits on the submillimeter flux. Our sample is shown in Table 1. 



Photometry observations were conducted using the SCUBA bolometer at the James 
Clerk Maxwell Telescope (JCMT) on Mauna Kea, Hawaii between February 2003 and Jan- 
uary 2004. The survey was only carried out in stable weather conditions when the precip- 
itable water vapor level was less than 2 mm. Zenith optical depths ranged from 0.1 to 0.2 at 
850 fim and 0.3 to 0.8 at 450 fim and the integration time spent on each source was between 
40 minutes to 3 hours. The resulting rms errors varied from 1 to 3 mJy beam~^ at 850 /xm 
and 10 to 30 mJy beam~^ at 450 //m. The pointing was checked via observations of bright 
quasars near each source after each major slew and the focus was adjusted every three hours 
on average, more frequently at times near sunrise and sunset. Calibration was performed 
by observations of Uranus, Neptune, and Mars and standard sources, CRL 618, CRL 2688, 
IRC+10216. Based on the agreement of the measured fluxes with those predicted for these 
calibrators, we estimate the photometry is accurate to within 10% at 850 /im and 30% at 
450 /xm. 

Additional CO (3-2) observations of HD 107146 and HD 104860 were made using the 
345 GHz receiver and 500 MHz acousto-optical spectrometer at the Caltech Submillimeter 
Observatory (CSO) in March 2003. The beamsize at this frequency is 15". The weather 
was very dry and stable resulting in system temperature that ranged from 450 K to 550 K. 
Spectra were taken in chopping mode with a throw of 120" at 0.321 Hz. The passbands 
were very flat over the frequency interval of interest and a single order baseline was removed 
from individual spectra before coadding and rebinning to 1 km s~^ channels. Pointing was 
checked via observations of IRC-l-10216 and was accurate to 2". Antenna temperatures were 
converted to the main beam scale using an efficiency ri„^i, = 0.65. No emission was detected 
in either source at an rms noise level of 5 mK (HD 107146) and 13 mK (HD 104860) per 
channel. 



2. 



Observations 



3. 



Results 



3.1. 



Submillimeter Continuum 



We detected 3 objects out of our sample of 13. The measured fluxes and upper limits 
for the entire sample are shown in Table 1. The errors and upper limits in this Table are 
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based on the rms noise only and do not include the calibration uncertainty discussed above. 
The results are also shown graphically in the spectral energy distributions (SEDs) plotted in 
Figures 1, 2, and 3. To construct the SEDs, we tabulated infrared fluxes from the 2MASS 
All-Sky Survey, the IRAS Faint Source Catalog (FSC), and additional ISO and submillimeter 
fluxes from the literature (Silverstone 2000; Carpenter et al. 2005; Liu et al. 2004; Habing 
et al. 2001). IRAS detections were color corrected based on the effective temperature of the 
star for 12 and 25 //m and (typically) for a 50 K disk at 60 and 100 /im(cf. HD104860, 
Figure 1, where 40 K was used). For those stars not in the FSC, upper limits are plotted 
based on the catalog completeness (0.2 Jy at 12, 25, 60 iim and 1.0 Jy at 100 //m, non-color 
corrected). The stellar photosphere is shown as a Kurucz model for the appropriate stellar 
type scaled to the 2MASS JHK fluxes. 

The brightest source detected at submillimeter wavelengths, HD 107146, is a G2 star 
with estimated age ~ 100 Myr. It was strong enough to map and found to be marginally 
resolved at 450 /xm. The maps and analysis of the SED were presented in Williams et al. 
(2004; hereafter Paper I). Briefly, the flt to the SED implied a dust temperature of 5 IK 
and a dust mass of O.IM®. Based on the measured properties of the source, we predicted 
that the disk would be unusually bright in scattered light. This was conflrmed in recent 
imaging of the system with the Hubble Space Telescope ( Ardila et al. 2004) . The other two 
submiUimeter detections in our survey are associated with F8 stars: HD 104860 with an 
estimated age of ~ 40 Myr and HD 8907 with an age ~ 200 Myr. 

For these two systems, we fit the excess above the predicted photosphere using a single 
temperature greybody with optical depth 

rA = l-exp[-(Ao/A>«] 

which has the desired asymptotic behavior, r = 1 for A ^ Aq and r = (Aq/A)'^ for A ^ Aq. 
The critical wavelength, Aq, was set to 100 fim for consistency with previous debris disk 
studies (Dent et al. 2000; Wyatt, Dent & Greaves 2003). Further, because the photospheric 
excess was significantly detected at only two wavelengths longward of 60 pLm in each case, 
we also fixed P — 1 based on the typical SED slope for debris disks (Dent et al. 2000) 
and protostellar disks (Beckwith et al. 2000). The excess was then fit by varying the dust 
temperature, Tdust, and scaling. The dust mass was determined from 

K'85oBi,{Tiiust) 

where Fs5o is the SCUBA 850 //m fiux, d is the distance to the source, ksso is the dust mass 
absorption coefficient, and Bi, is the Planck function. 

The value of /3 ~ 1 for disks is less than that observed in the ISM (/3 ~ 2) and can 
be produced by a population of large grains (e.g., Miyake & Nakagawa 1993; Pollack et al. 
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1994). Independent evidence for large grains in debris disks comes from the lack of spectral 
features at mid-infrared wavelengths indicating that the emission arises from grains > 10 /xm 
in size (Jura et al. 2004) and is consistent with a coUisional origin for the dust (e.g., Krivov 
et al. 2000). 

Pollack et al. (1994) show that for spherical grains at 100 K, < 1 is only achieved for 
radii > 1 cm. The resulting values of k,s50 range from 0.32 to 0.62 cm^ g~^. These are a factor 
of 3 — 5 less than the canonical value, k, = 1.7 cm^ g^^, that is commonly used in studies 
of debris disks dating back to Zuckcrman & Becklin (1993). On the other hand, porous, or 
fractal, dust grains can have cross-sections that are more than an order of magnitude higher 
than spheres per unit volume (Wright 1987). Due to the uncertainty in the dust composition 
and structure, we also adopt k — 1.7 cm^ g~^ for ease of comparison with previous results 
but we note that the masses may be 3 — 5 times higher than the values assumed here. 

HD 104860 was undetected by IRAS at mid- to far-infrared wavelengths, but we mea- 
sured emission in excess of the photosphere at both 850 /im and 450 /im indicating a cold, 
massive disk. The SED is shown in Figure 1 where the photospheric excess is fit by a range of 
/? = 1 greybodies with temperatures between 19 and 42 K. The best fit has T^ust = 33 K and 
Mciust = 0.16 M^. Table 2 summarizes the corresponding range of dust masses and excess 
luminosities. The lack of an IRAS detection and the low temperature might suggest that 
the submillimeter flux is due to Galactic cirrus emission but the high latitude of the source, 
b — 50°, makes this unlikely. The low temperature of the dust surrounding HD 104860 is 
similar to that of three other cold disks that have been reported around companions to early 
type stars (Wyatt et al. 2003). 

HD 8907 was detected by IRAS at both 12 ^m and 60 /xm and by ISO at 60 /xm and 
87 yum (Silverstone 2000). The 12//m flux is photospheric, whereas the longer wavelength 
fluxes represent a strong excess. In the submillimeter, an excess is detected at 850 /im and 
marginally {2a) at 450 iim. The P — 1 greybody flt to the excess above the photosphere 
is tightly constrained by these points with a temperature T^ust = 48 K, and mass M^ust — 
0.036 Me (Figure 2). 

The SEDs of the remaining ten non-detections are shown in Figure 3. Three objects, 
HD 17925, HD 35850, and possibly HD 166435 show mid-infrared excesses above the stellar 
photospheres, but lack detections longward of 100 fim. The 3-a upper limits to the 850 
flux and inferred disk masses are presented in Table 1. The flux upper limits were converted 
to mass upper limits using an average temperature for debris disks around solar type stars. 
Since stars with luminosities < 10 Lq have dust temperatures within a narrow range 33 
to 56 K, averaging 50 K (Wyatt et al. 2003; Sheret et al. 2003; Sylvester et al. 2001; Liu et al. 
2004; Greaves et al. 2004b), we therefore adopted an average temperature of {Tdustj = 50 K 
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in estimating the mass upper limts. Note that all the upper limits are less than or equal to 
the most massive disk in our sample, HD 104860, and five are less than the least massive 
disk, HD 8907, and therefore these non-detections provide statistically useful information. 
The range in source distances of course enters into the mass upper limits, but the sensitivity 
of the observations is nevertheless adequate to show that there is a range of disk masses 
present around stars of similar age. For example, for HD 17925 at a distance of 10 pc and 
an age of 100 Myr, our mass upper limit is O.OOSM^. In comparison, the more distant 
source HD 107146 (at 28.5pc) has a similar age and a mass that is 20 times larger. These 
results provide further evidence for intrinsic diversity in the evolution of the circumstellar 
dust content of solar-type stars (cf. Meyer et al. 2004). 

3.2. CO (3-2) 

We can convert the upper limits on the CO(3-2) integrated intensities for HD107146 and 
HD 104860 to upper hmits on the CO gas mass assuming LTE and an excitation temperature 
for the gas (e.g. Scovillc ct al. 1986). For HD107146. wc assume an excitation temperature 
of Tex = 50 K based on the estimated dust temperature for the system (Paper 1) and a 
FWHM for the line of 1.3 — 2.3 km s^^, which is appropriate for gas at ~ 30 — 100 AU 
in orbit around a ~ IM0 star at an inclination of 25° (Ardila et al. 2004). For the more 
conservative limit corresponding to a FWHM of 2.3 km s~^, the 3a upper limit to the mass 
of gaseous CO is Mco < 2 x 10"^ M®. For HD104860, we assume T^^ = 30 K based on 
the estimated dust temperature for the system (Table 2). Since the low dust temperature 
is consistent with dust emission from large disk radii (> 100 AU), 2 km s"~^ FWHM is a 
reasonable line width. With these assumptions, the 3a upper limit to the mass of gaseous 
CO is Mco < 1-3 X 10"^ Me. These hmits vary by less than 50% for excitation temperatures 
in the range T^^ = 25 - 100 K. 

A critical issue is the conversion from CO to total gas mass. Processes that can reduce 
the abundance of gas phase CO include condensation on grains and photodissociation (e.g.. 
Dent et al. 1995; Kamp & Bertoldi 2000). If the upper limit on the CO mass is to be 
consistent with a primordial gas-to-dust ratio of 100 (i.e., a total gas mass of lOM®), these 
processes must be extremely efficient given the large dust mass and our upper limit on the CO 
mass. For example, for the dust mass of O.IOM^ measured for HD107146 (Paper 1; see §3.1), 
a gas-to-dust ratio of 100 implies a CO to total hydrogen abundance of nco/n-a^ = 2 x 10~^ 
(assuming a 10% fraction of Helium by number), several orders of magnitude below the 
interstellar value of riHa/'T'CO — 10^ (e.g., van Dishoeck et al. 1993). 

For the condensation of CO on grains, an important issue is the substrate onto which 
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CO is adsorbed. The binding energy of CO onto water ice is 1740 K, whereas the CO-CO 
surface binding energy is significantly lower, 960 K (Sandford & Allamandola 1988, 1990). In 
their study of well-known debris disks, Kamp & Bertoldi (2000) assumed CO adsorption onto 
a water ice substrate, in which case where the critical grain temperature for condensation 
turns out to be ~ 50 K. In contrast, studies of CO condensation in protoplanetary disks 
have typically assumed a CO substrate where the grain condensation temperature is ~ 20 K 
(e.g., Aikawa et al. 1996). The relevant substrate likely depends on the thermal history of 
the grains in the disk. 

For example, if the grains are imagined to be heated and rapidly cooled, gas phase water 
and CO will condense together, with CO binding to water ice. In contrast, if the grains are 
heated and slowly cooled, water will condense first, followed by CO. In the latter case, the 
abundance of CO in the gas phase relative to the number of grains is sufficient to establish 
a CO substrate. That is, for a gas-rich disk in which the gas-to-dust mass ratio is 100, with 
all of the dust in grains of radius a, and nco/^H ~ 10~^, there are ~ 10^'^ CO molecules 
per 10/um dust grain. In comparison, there are Aira'^ns binding sites available on the grain, 
where ~ 1.5 x 10^^ cm~^ (Tielens & Allamandola 1987), or ~ 2 x lO^'' sites per monolayer. 
Thus, the condensation of only a very small fraction (< 1%) of the CO is needed to establish 
a CO substrate. In this case, because of the low CO-CO surface binding energy, the critical 
temperature for condensation will be lower than that estimated by Kamp & Bertoldi (2000). 

For HD107146, the measured temperature of the emitting grains is high (~ 50 K), 
comparable to the evaporation temperature for even CO on water ice. If the emitting grains 
are large and therefore at the blackbody temperature, condensation of CO on grains is not 
significant. Such grains would be located close to the star, at ~ 30 AU. However, if the ~ 50 K 
grain temperature characterizes the emission from a population of small (non-blackbody) 
grains located farther from the star, a larger, cooler population of grains, if also present, 
could potentially deplete CO from the gas phase. A distance of 100 AU is a reasonable 
maximum distance for the emitting grains given the spatially resolved continuum emission 
at 450/im (Paper I). At this distance from the star, grains > lOO/xm will be at the blackbody 
temperature of 28 K. Significant condensation will begin to occur when the adsorption rate 
onto grains exceeds the evaporation rate. For condensation onto a CO substrate, evaporation 
will keep CO from condensing further up to a gas phase CO density of nco < 3 x lO^cm"^. 
This density limit exceeds not only the density expected at the midplane of the minimum 
mass solar nebula, but also the much lower densities expected in the lower mass HD 107146 
disk. Thus, a gas-rich disk is not expected to experience significant CO condensation for 
grains located anywhere from the inner radius of ~ 30 AU out to ~ 100 AU. 

For HD104860, the temperature of the emitting grains is lower, ~ 30 K, and the uncer- 
tainty in temperature allows the possibility of a temperature as low as 19 K (Table 2). Given 
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the possibility that the emission arises from a population of relatively small (non-blackbody) 
grains and that larger, cooler grains may be present, it is difficult to rule out the possibility 
of significant CO condensation in this case. 

Another process that can reduce the gas phase CO abundance is photodissociation. 
Kamp & Bertoldi (2000) have shown that FUV irradiation by the star and the interstellar 
radiation field can produce significant underabundances of CO in the disks around Vega 
and 13 Pic. Using these results, they argue that the non-detection of CO in these systems 
may be consistent with a primordial gas-to-dust ratio in these systems. Other studies of 
debris disks surrounding both A-type and solar-type stars have previously argued that CO 
has been significantly dissociated and that the upper limits on CO mass are consistent with 
primordial gas-to-dust ratios (Dent et al. 1995; Greaves et al. 2000). A later spectral-type 
star such as HD107146 is beheved to be a negligible source of dissociating fiux compared to 
the interstellar radiation field (Kamp & Bertoldi 2000; Greaves et al. 2000). 

Following Dent et al. (1995) and Greaves et al. (2000), we can use the photodissociation 
models of van Dishoeck & Black (1988) to estimate the extent to which CO will be dissociated 
in the HD 107146 disk. These models show that the relative abundance of CO is a function 
of both the density within and the column density though the medium. For a gas-to-dust 
ratio of 100, the expected disk gas mass is lOM^. This corresponds to a vertical hydrogen 
column density of ~ 4 x 10^^ — 4 x 10^^ cm~^ for gas uniformly distributed in a disk 30-100 
AU in radius, and a mean density of nn ~ 10'' — 10^ cm~^ assuming a thermal scale height. 
At such large column densities, CO molecules are expected to be well-shielded. The models 
of van Dishoeck & Black (1988) predict that at these column densities the ratio of the CO 
and hydrogen column densities is > 10~^, several orders of magnitude larger than the ratio 
of 2 X 10~^ required for a primordial gas-to-dust ratio. 

Thus, our CO mass limit appears to completely rule out a primordial gas-to-dust ratio 
for HD107146. How low might the actual gas content be? Once the total hydrogen column 
density of the disk begins to drop below lO^^cm"^, the CO molecules are exposed to a 
higher photodissociating fiux and the CO abundance relative to hydrogen declines below the 
well-shielded value. For a gas-to-dust ratio of order ~ 1, we would be in the regime where 
the expected hydrogen column density and hydrogen number density are A^h ~ 10^'' cm^^ 
and nn ~ 10"'' — lO'^cm"^, respectively. Model T6C of van Dishoeck & Black (1988; nn ~ 
IC^ cm^'^) provides the closest match to the expected range of density. At hydrogen column 
densities A^h ~ lO^'^cm"^, the model predicts that the ratio of CO and H2 column densities 
is A^co/A^Ha ~ 10~^ or a mass ratio of Mco/Mn^ ~ 10"^. Thus the model implies that our 
upper limit on the CO mass of 2 x 10~^M^ corresponds to an upper limit on the total gas 
mass of ~ O.O2M0, or a gas-to-dust ratio of < 1, significantly below the primordial value. 
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One of the possible flaws in this argument is that the van Dishoeck & Black models used 
to interpret the results were developed to address the conditions in interstellar clouds rather 
than those in circumstellar disks. So, for example, they describe lower density conditions 
than may apply to HD 107146. They also do not include the possible effects of additional 
processes, such as advanced grain growth and stellar X-ray and UV excess irradiation, on 
the chemical abundance structure of the medium. We have also assumed that the gas and 
dust temperatures are equal, an assumption that should be examined theoretically. An 
improved analysis would be possible with the recent models of Gorti & HoUenbach (2004) 
which include such processes and more directly address the thermal and chemical structure 
of residual gas disks. This is a topic for future investigation. 



4. Discussion 

4.1. Dust Disk Radii and Masses 

Single temperature greybody fits like those employed here are commonly used to inter- 
pret sparsely sampled SEDs of debris disk sources. Taken at face value, such a fit implies that 
the emission arises from a limited range of disk radii, in particular that the dust distribution 
has an inner hole. An inner hole in the dust distribution could arise in several different ways. 
For example, in a disk composed of icy grains, grain destruction at the water ice evaporation 
radius can produce an inner hole in the dust distribution. While this situation may apply 
to disks surrounding luminous (A- and B-type) stars that have dust temperatures ~ 100 K, 
similar to the water ice sublimation temperature (e.g., HR4796A — Jura et al. 1998), the low 
dust temperatures of the sources detected in our survey (30-50 K) suggest that this situation 
is not relevant here. 

Primordial gas-rich disks may also develop optically thin regions as a result of grain 
growth and planetesimal formation. Since collisional evolution is more rapid at small radii, 
optically thin inner holes arc expected to develop. In the gaseous outer disk, large dust 
masses can be maintained since gas drag can prevent grains from developing eccentric orbits 
and substantial relative velocities, and destructive collisions are thereby minimized. Such 
gas-rich outer disks with optically thin inner holes have been proposed to explain the large 
excesses detected in young disks systems such as TW Hya (~ 10 Myr; Calvet et al. 2002) 
and CoKu Tau/4 (~ 1 Myr; D'Alessio et al. 2004). Dusty disks that possess a residual gas 
component {Mg^s ^ lOMdust) can also develop ring- like structures as a result of gas-grain 
coupling in the presence of stellar radiation pressure. This situation has been proposed 
to explain the ring-like structures seen in systems such as HR4796A and HD141569 (e.g., 
Takeuchi & Artymowicz 2001). 
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The theoretical expectation that gaseous outer disks photoevaporate on short timescales 
(~ 10 Myr; HoUenbach, Yorke, & Johnstone 2000) and the lack of observations to the 
contrary both suggest that the above situations are unlikely to apply to our sample (age = 
10-200 Myr), although they arc difficult to rule out. For solar mass stars, irradiation by the 
central star is thought to be capable of photoevaporating away the gaseous component of 
the outer region (r > 10 AU) of a minimum mass solar nebula on a timescale of 10 Myr, if 
the Lyman continuum flux from the star remains at classical T Tauri levels for this length of 
time (HoUenbach et al. 2000). Since the Lyman continuum flux level is unlikely to continue 
at such levels for longer than 10 Myr, it may be difficult to dissipate more massive disks 
through photoevaporation by the central star alone. TW Hya may be one such example of 
a ~ 10 Myr old star with a substantial outer gas disk that has survived an earlier phase 
of energetic photoevaporation (Kastncr ct al. 1997; Qi ct al. 2004 and references therein). 
Other processes such as viscous dissipation may then dominate the dissipation of the disk. 

Observational support for the rapid dissipation of gaseous outer disks comes from the 
results of searches for gas in 1-10 Myr sources (e.g., Zuckerman, Kastner, & Forveille 1995) 
and known debris disk systems (e.g.. Greaves et al. 2000) The limited data available for our 
sample (sec. 3.2) indicate an upper limit to the gas mass for HD107146 that is less than the 
dust mass in the system, but there is no significant constraint on the gas content of the rest 
of the sample. Even for HD107146, although our results apparently rule out a primordial 
gas-to-dust ratio for the system, the upper limit on the gas mass may still allow gas-grain 
coupling and radial migration of grains to play a role in sculpting the dust disk. 

The possibility that these studies, which primarily do not trace H2 directly, have un- 
derestimated the disk gas content has been raised by the surprising detection with ISO of 
abundant (~ IMj) H2 in ~ 20 Myr old debris disk systems (Thi et al. 2001). Spitzer offers 
an opportunity to confirm these measurements and to make similar measurements for a wider 
range of debris disk sources. Some uncertainty already surrounds the ISO results given since 
ground-based searches for mid-IR H2 emission from younger T Tauri stars do not confirm 
the ISO detections (Richter et al. 2002; Sheret et al. 2003; Sako et al. 2005). In addition, 
recent Spitzer observations do not confirm the reported H2 line fiuxes for jS Pic (C. Chen, 
personal communication). Whether the imphed lower gas masses are representative of the 
general situation for debris disk systems at 20 Myr and what constraints are placed by the 
observations on the disk gas mass are interesting issues for the future. 

Given the existing uncertainties regarding the gas dissipation timescale in disks, it is 
difficult to rule out the possibility that sources in our sample may have significant gaseous 
outer disks; gas-grain interactions may therefore contribute to the structure that is deduced 
from the observed SED. For the rest of the paper we assume, consistent with convention, 
that at the relatively advanced age of the sources in our sample (10-200 Myr) httle residual 
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gas remains in the outer disk region. Sensitive future searches for gas in these or similar 
systems, such as those to be made with Spitzer as part of FEPS, may spark a reexamination 
of this issue. 

In the absence of a substantial gaseous component to the disk, grain lifetimes are affected 
by both Poynting-Robcrtson (PR) drag and particle collisions. Collisions are expected to 
dominate the destruction of grains at high disk masses (e.g., Dominik & Decin 2003; Wyatt 
2005). As described by Backman & Paresce (1993), the grain mutual collision time is tcoU ~ 
P/8a where P is the orbital period and a is the fractional surface density (the surface 
filling factor) in grains. This estimate assumes that grains are on orbits that are sufficiently 
inclined that the grain encounters the full surface density of the disk approximately twice per 
orbit as it oscillates through the disk plane. The estimate also includes tangential collisions 
as destructive events, i.e., the coUisional crosssection is 4 times the crosssectional area of 
individual grains. For a total dust mass M^usti spread out over an area TrR'^^g^, internal grain 
density pi, and average grain radius a, 

tcou = 3000yr (i?d„,t/30AU)'/' {Mdust/O.lM^y' (a/lmm) . 

Thus, for the sources in our sample with detected dust disks {M^ust ~ O.IM®; Table 1), 
the collision timescale is significantly shorter than the age of the systems, and the observed 
grains must be replenished by collisions between larger bodies. 

In comparison, the PR drag time for blackbody grains is 

tpR = 640 Myr (a/lmm) {Raust/^^AUf {LJLq) 

where is the stellar luminosity. For the millimeter-sized grains deduced from the submil- 
limeter observations, collisions therefore dominate as the destruction process down to dust 
masses Mdust < IO^^Mq, well below the dust masses inferred for the systems studied here. 
At dust masses above ~ lO-^M®, grains are destroved by collisions before PR drag can cause 
them to migrate very far from where where they were created. As a result, the region of the 
disk that contains the coUisional debris is also the region of the disk in which the colliding 
planetesimals that produce the debris are located. 

We can estimate the orbital radii of the coUisional debris for the sources detected in 
our sample assuming that the grains responsible for the emission are large and absorb and 
emit as blackbodies. The ~ 50 K dust temperature of HD107146 and HD8907 translates 
into grain orbital distances of approximately 30 and 48 AU respectively. The lower (~ 30 
K) temperature of the dust in the HD104860 system locates the coUisional debris at ~ 110 
AU. We therefore infer the presence of a population of colliding planetesimals at these radii. 
The presence of a large planetesimal population in these young solar analogue systems at 
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orbital radii similar to that of the Kuiper Belt in our solar system (e.g., Jewitt & Luu 2000, 
PPIV) suggests their utility in understanding the early evolution of Kuiper Belt systems. 

4.2. Comparison With Models 

Several recent models for the coUisional evolution of planetesimal disks suggest possible 
interpretations for these data. Dominik & Decin (2003) describe the evolution of a ring 
of planetesimals located at a fixed distance from the star. Assuming that the ring is in 
equilibrium, in the sense that grain production and destruction rates balance, the mass in 
small grains evolves as Mdust oc if grain destruction is dominated by mutual collisions 
and Mdust oc if grain destruction is governed by PR drag (see also Spangler et al. 2001). 
The former case would apply at early times when the disk mass is large and collisions are 
frequent; the latter would dominate at late times when collisions become infrequent enough 
the PR drag begins to dominate. The model is flexible in that it does not specify how the 
eccentricity and inclination distributions of the planetesimals are determined, i.e., how and 
when the population is stirred. 

In a complementary, more prescriptive model, Kenyon & Bromley (2004) follow the evo- 
lution of a disk of planetesimals undergoing planet formation. They find that once embedded 
objects grow to > 1000 km in size, they stir up the surrounding planetesimals, generating a 
collisional cascade and accompanying debris. That is, the planetesimal disk is "self-stirred" 
by embedded protoplanets. Since the evolution of the disk is dominated by collisions, the 
growth timescale scales as the collision time tcou oc P/S where P is the orbital period and 
S is the disk surface density in planetesimals. Since both the orbital speed and disk surface 
density decline with disk radius, the growth time is an increasing function of disk radius. 
Thus, a wave of planet formation propagates outward in the disk: as time progresses, the 
formation of embedded protoplanets excites collisions among the less massive bodies, thereby 
generating dust debris at successively larger characteristic radii. 

In this model, the total mass in small grains in the disk remains remarkably constant 
with time (for perhaps as long as 1 Gyr), and subsequently decays as oc where n — 1 — 2. 
In contrast, the reprocessed luminosity emitted by the collisional debris begins to decline at 
a much earlier time (after only 1-10 Myr). This is because as the wave of planet formation 
moves outward, grains of a given size subtend increasingly smaller solid angles as they are 
located at larger distances from the star. Although the models calculated by Kenyon & 
Bromley are for more massive stars (M* = SM©) than those in our sample, we can use 
the scaling of tcoii oc P/E to estimate the behavior that would be obtained for lower mass 
systems. 
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Comparing our results with these models, we find that the deduced size of the dust disks 
for the detected systems and their nominal stellar ages are roughly consistent with predictions 
for the evolution of dust production in collisional planetestimal disks. The Kenyon & Bromley 
(2004) model predicts that in a disk with a mass comparable to that of the minimum mass 
solar nebula that surrounds a solar mass star, it takes ~ 300 Myr for planet formation to 
reach 40 AU and to produce a local maximum in the dust production rate at that distance. 
This is comparable to the sizes (30-50 AU for HD8907 and HD 107146) and ages (200 and 
100 Myr respectively) of the sources studied here. 

The derived dust masses for the sources detected here are 0.03 — O.IGM®. These dust 
masses are 2-8 times the mass in small grains that is predicted by Kenyon & Bromley (2004) 
to reside in a disk with the mass of the minimum mass solar nebula that is undergoing planet 
formation (0.01 — 0.02M®). In the Kenyon & Bromley (2004) study, the mass in collisional 
debris was found to be remarkably independent of the bulk properties of the planetesimals 
and initial conditions other than disk mass. Instead, the dust mass was found to depend 
primarily on (and to be proportional to) the initial mass of solids in the disk. Interpreted in 
this context, the larger disk masses measured here suggest that they arise from more massive 
versions of the minimum mass solar nebula. 

The Kenyon & Bromley (2004) calculation provides a way in which to deduce the mass 
of planetesimals that currently resides in the disk (see also Dominik & Decin 2003). Their 
results predict that at an age of 100 Myr ~ 4 x 10"*^ of the mass in solids is in the form 
of grains 1 mm and smaller. Thus the measured mass in small grains implies a current 
mass in parent bodies of ~ SSM® for HD8907, ~ 250Me for HD107146, and ~ 400Me for 
HD104860. They also find that over the age range of the sources in our sample, the mass 
in small dust grains is roughly half that of the initial mass of sohds in the disk. Thus the 
initial mass in solids in the outer planetesimal disks of the detected sources would be twice 
as large (~ 170 — SOOAfe). For comparison, the mass of the Kuiper Belt in our solar system 
today is a mere ~ O.lA'/0; in primordial times, the minimum mass solar nebula in the entire 
radial range 30-150 AU contained ~ lOOM® and several 10s are thought to have been 
needed in the 30-50 AU range in order to form Pluto and QBi-type 100m KBOs (Stern & 
Colwell 1997; Kenyon & Luu 1998). 

Since our sample focuses on solar mass stars, it is interesting to ask what relevance these 
stars might have to our solar system. The large dust masses measured for the submillimeter 
sources detected in our study indicate initial disk planetesimal masses 2-10 times that of the 
minimum mass solar nebula. Such massive disks favor the production of giant planets like 
those in the solar system. Massive disks make it easier to form Jupiter with the possibly 
short hfetime of the gaseous disk (Thommes et al. 2003) . So these systems may be reasonable 
analogues of the young solar system. 
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The inferred masses of the planetesimal disks in the current epoch are massive enough 
to play a role in determining the architectures of planetary systems. Planetesimal disks ~ 
lOOAf® are expected to induce giant planets to migrate significantly on ~ 100 Myr timescales 
(Hahn & Malhotra 1999), since the clearing of a residual planetesimal disk by planets involves 
a significant exchange of orbital energy and angular momentum. In the case of our solar 
system, the interaction between a 10— lOOM® residual planetesimal disk and the giant planets 
is believed to have caused Jupiter to migrate inward and Saturn, Uranus and Neptune 
to migrate outward (e.g., Hahn & Malhotra 1999). In particular, the significant outward 
migration of Neptune (by ~ 7 AU) is thought to have sculpted the Kuiper Belt and pushed 
it out to the orbital distance range that it currently occupies (40-50 AU; Hahn & Malhotra 
1999; Levison & Morbidelli 2003). The disks detected here appear to be sufficiently massive 
to induce similar migration among any giant planets that may be present. 

When interpreted (literally) in the context of the Kenyon &; Bromley (2004) model, the 
coUisional debris that we have observed surrounding the ~ 100 Myr sources in our sample 
suggests that we are witnessing the formation of ~ 1000 km bodies at distances of 30-70 
AU around nearby stars. For HD 107146 and HD8907, which have debris disks at ~ 30 — 50 
AU, the large bodies that are forming may correspond to proto-Plutos. These disks are 
comparable to the size of the Kuiper Belt in the solar system, which appears to have an edge 
at ~ 50 AU (Allen, Bernstein, & Malhotra 2001; Chiang & Brown 1999). The possibihty of 
a primordial Kuiper Belt extending beyond 50 AU has been discussed in the literature (e.g.. 
Stern 1995), and indeed HD104860 is one example of a young (40 Myr) solar-mass star with 
a planetesimal disk that appears to extend much beyond known extent of the Kuiper Belt 
in our solar system (to Rdust ~ 110 AU) and is massive enough at that distance to enable 
the formation of a massive perturbing object. 



4.3. The properties of disks detected at 850 /xm 

Since the inferred properties of the debris disks (planetesimal and initial disk masses) 
and potential implications discussed above rely on the calculations of Kenyon & Bromley 
(2004), one might wonder about the applicability of these models to the sources that we 
observed. For example, a possible limitation of the Kenyon & Bromley study is that it does 
not include the possible formation of giant planets at small radii and their impact on the 
coUisional evolutionary history of an outer planetesimal disk. The presence of shepherding 
giant planets has often been invoked to explain the ring-like structures, observed or inferred, 
in disk systems over a range of ages, from ~ 1 Myr (e.g., CoKu Tau 4; D'Alessio et al. 2004) 
to ^ 100 Myr (e.g., e Eri; Greaves et al. 1998). In particular the likelihood of a diversity of 
giant planet architectures, similar to that seen in the extrasolar planet population (Marcy 
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& Butler 1998), may lead to a diversity of coUisional evolutionary histories that depart from 
those predicted by self-stirred models. One way to explore the extent of such resultant 
diversity and the applicability of the models to observations of outer planctcsimal disks is 
to compare the model predictions with the collective properties of debris disks detected at 
submillimeter wavelengths. 

We therefore collated the submillimeter excess sources from the literature and combined 
them with the results of our survey. In practice, almost all detections of debris disks at 
sub-millimeter wavelengths have been made with the SCUBA camera on the JCMT at 
850 /xm(Greaves et al. 2004a, 2004b; Wyatt et al. 2003, 2004; Sheret et al. 2004; Liu et al. 
2004). We henceforth refer to this sample as the 850 /xm disks. In the following sections we 
examine compare the masses, radii, and ages for the sample with the predictions of recent 
models for the coUisional evolution of planetesimal disks. 

4-3.1. Disk Mass Evolution 

Figure 4 plots the dust masses and ages of submiUimeter excess sources from the litera- 
ture combined with the detections reported here. For these sources, the SCUBA fluxes were 
converted to masses assuming k — 1.7 cm^ g~^ and dust temperatures from the literature. 
Non-detections are indicated in Figure 4 as 3-0" upper limits; these are calculated assuming 
a dust temperature of 50 K. The sources shown cover a range of spectral types, from B7 to 
M4. The resulting plot is similar to those that have appeared previously in the literature 
(e.g., Wyatt et al. 2003; Liu ct al. 2004). Our detections follow the trend established by the 
earlier measurements, which shows a decline in the dust mass as a function of age. Such a 
decline is expected to result from the coUisional erosion of a planetesimal disk. 

Several significant selection effects complicate the interpretation of this plot. For exam- 
ple, the detected sources are mostly early type stars, which is likely a result of sensitivity 
considerations. If these more massive stars form with more massive disks, they will have 
preferentially larger dust masses at any given age, making them easier to detect. The higher 
luminosity of the central star also enhances the detectability of any surrounding dust. An- 
other selection effect is that there are few sources < 10 Myr in age that are near enough to 
be readily detectable at low masses. (The sources in the TW Hya and (3 Pic moving groups 
are notable exceptions.) Since such sources are typically located at larger distances than the 
older sources that have been detected, the lack of sources in the lower left region of the plot 
may simply result from sensitivity limitations. Thus, although there are no sources with low 
measured dust masses at 10 Myr age, there are many upper limits. 

The upper right region of the plot is also not well populated, a property that seems 
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difficult to attribute to a selection effect. In addition, some of measurements made here 
(in the 10-100 Myr range) faU below the mean mass-age relation. This suggests that the 
detected sources represent the upper envelope of a broader distribution. This is consistent 
with the interpretation discussed earlier that the disks detected here arise from more massive 
versions of the minimum mass solar nebula. Much more massive disks are not expected since 
they would be gravitationally unstable. 

To investigate the constraints that these data place on the evolution of debris disks, in 
particular the information contained in the many non-dctcctions, we carried out a survival 
analysis (Feigelson & Nelson 1985) using the ASURV Rev 1.2 package (LaValley, Isobe, & 
Feigelson 1992). We divided the sample into a 'young' and 'old' group based on a dividing 
age, to, and then computed the probability that they were drawn from the same population. 
The probabilities of similarity decreased from 7 - 12% at to = 100 Myr to 2 - 3% at 200 Myr 
and down to 1 — 2%^ at 500 Myr. Thus, there appears to be a significant difference between 
the 'young' and 'old' groups at a dividing age of to ~ 200 Myr. That is, the mass in small 
grains decreases significantly on a ~ 200 Myr timescale. 

This result is complementary to that obtained by Carpenter et al. (2004). Using survival 
analysis, they compared the submillimeter excess detections and upper limits obtained for 
a younger sample of FEPS sources with the submillimeter properties of sources in Taurus. 
They found no statistically significant difference between the Taurus disk population and 3- 
10 Myr systems, somewhat weaker evidence for a difference with the 10-30 Myr population, 
and a clear difference between the Taurus disk population and 30-100 Myr systems. The 
results presented here suggest the decline to lower dust masses continues at ages beyond 
~ 200 Myr. 

Although the small size (and heterogeneous nature) of the submillimeter sample makes 
it difficult to constrain the evolutionary behavior of the sample in greater detail, the dashed 
hne in Figure 5 shows that the measurements can be fit by an overall 1 / (Age) dependence 
of dust mass vs. age with a significant spread about this trend. This overall trend is the 
behavior expected for a ring of planetesimals in which grain destruction is governed by 
collisions (Dominik & Decin 2003). This slope is also consistent with that predicted in the 
Kenyon & Bromley (2004) models for ages > 10 Myr. Some of the spread in the observed 
dust mass at a given age may be due to the range of stellar spectral types present in the 
sample. For example, Dominik & Decin (2003) show that a range of spectral types AOV to 
K3V introduces a spread of about an order of magnitude in fractional dust luminosity at 



The precise probabilities depended on the particular statistical test used and the weighting that they 
give to the upper limits. 
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an age of 1 Gyr, a range comparable to the spread found here. Additional heterogeneity in 
the sample, such as a range of initial disk masses or a diversity of giant planet formation 
histories, would also contribute to a range in disk mass at any given age. 

A caveat to these results is the possibly significant uncertainty in the estimated ages 
for the stars in the sample. The ages in Table 1, which derive from those estimated by 
FEPS, arc typically uncertain by factors of 2-3. Similar errors are likely to characterize 
the ages of the sources from the literature. Note, however, that larger age uncertainties 
may characterize some sources. For example, subsequent to their selection for this study, 
HD88638 and HD166435 were both dropped from the FEPS sample because their Li strengths 
indicated ages older than the ~ 10 Myr age estimated initially for these sources. So while 
typical age errors of a factor of 2-3 would not significantly affect our conclusions, this issue 
may need to reexamined if larger errors are eventually found to be common. 

Since the survival analysis formalism used above does not account for errors in the vari- 
ables, we tested the sensitivity of the inferred decline in mass with age by varying individual 
stellar ages by up to a factor of 3. We found that only when several systems crossed over 
the 200 Myr timeline did the statistics change substantially. The effect was to broaden the 
dividing age to into an age range of 70 — 300 Myr over which an "old" population is signif- 
icantly less massive than a "young" population. The two populations and the approximate 
border between them are shown in the two-tone greyscale in Figure 4. 



4-3.2. Disk Temperatures and Inner Radii 

The Kenyon & Bromley (2004) models make the simple prediction that the characteristic 
dust radius increases with age in a given system if self-stirring dominates the collisional 
evolution of planetesimal disks. We therefore compared the stellar ages with the dust radii 
inferred for the sample. The dust temperatures of the sources, either as measured in section 
3 or from the literature, were converted to characteristic dust radii assuming blackbody 
grains: Rdust,bb = 7.8 x 10^ A\J (L/ L^f-^ 

^dusf Stellar ages were taken from the literature. 

Figure 5 plots the derived disk radius as a function of stellar age. No obvious correlation 
is found for the sample as a whole. Since the large range in stellar spectral types may obscure 
a trend, it would be preferable to compare the properties of a sample spanning a smaller 
spectral type range. Among the A stars in the sample no strong trend is found although 
sample statistics arc admittedly poor. We might attempt to include all of the sample in a 
comparison by making an assumption about how the dust evolution depends on the stellar 
mass. If the time for debris production to reach a given radius Rd is related to the collision 
timescale tcoii oc P/E where the orbital period P oc R^J'^M^^^'^ and the surface density in 
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planetesimals follows the slope inferred for the solar nebula T, oc , then tcoii oc -R^M* ' . 
Assuming tage ~ tcoii, we would expect tage oc i?^M*~^^^. Since the stellar mass of the sample 

1/2 

varies by only 6 in the sample, scaling by produces only a minor variation on Figure 5. 
Another hypothesis might be that the initial planetesimal disk mass scales with stellar mass 
so that S oc M^R^ .In this case tcou ~ R^M^ . The stronger dependence on M* tightens 
up the distribution slightly but the spread is still large and no strong trend is apparent. 

There arc several caveats to these results. For one, we have assumed that the grains 
responsible for submillimctcr excesses arc large enough that they emit as blackbodies. This 
assumption is supported by recent results from Spitzer which show that debris disks typically 
do not show emission features, suggesting that the grains responsible for the excess are large 
> 10//m (Jura et al. 2004). Similar spectroscopic results for the specific objects studied here 
would put our assumption on firmer footing. 

Perhaps more significantly, sparsely sampled SEDs do not allow detailed fits that inde- 
pendently constrain the dust temperature and (3. (The spectral slope /9 is a function of grain 
composition and shape as well as the grain size distribution.) As a result, we have typically 
assumed j3 — 1 m. our fits to the dust temperatures. This may be an important source of 
error. For example, for HD 107146, where the SED is particularly well constrained, the best 
fit SED implies T = 51 K and /3 = 0.7 (Williams et al. 2004), compared to a temperature of 
43 K that is obtained assuming j3 = 1. This ~ 20% difference in temperature translates into 
a factor of 40% error in dust radius. This is a significant difference although much smaller 
than the dispersion in Figure 5. 

Finally, one might ask how the inferred dust radii compare with the spatially resolved 
sizes of debris disks where available. The radial extent of the excess emission, when spatially 
resolved at long wavelengths, generally agrees within a factor of 1.5 with the dust radii es- 
timated from the SED. Examples of such systems include Fomalhaut (Holland et al. 1998), 
HR4796 (Jayawardhana et al. 1998; Koerner et al. 1998), and Vega (Wilner et al. 2002). 
In contrast, the sizes of debris disks observed at short wavelengths (optical through near- 
infrared) are typically larger than the sizes inferred for the same systems based either on 
their SEDs at longer wavelengths or through resolved millimeter or submillimeter imaging. 
Examples of such systems include Vega (Su et al. 2005), j3 Pic (Holland et al. 1998; Kalas & 
Jewitt 1995) and AU Mic (Liu et al. 2004; Kalas et al. 2004). In such cases, the short wave- 
length scattered light probably traces a small grain population that is either being ejected or 
placed into eccentric orbits through a combination of radiation pressure, corpuscular (stellar 
wind) drag, and dynamical scattering with massive planets (e.g., Moro- Martin & Malhotra 
2002, 2003b). 

Thus, the available measurements suggest that the dust radii determined from the SED 



-19- 



assuming blackbody grains reasonably estimates the radius of the emitting grains. This could 
be confirmed by spatially resolving the submillimeter excess emission in a larger number of 
sources. When such measurements are available it will be interesting to examine the true 
behavior of dust radius as a function of the age of the system. At the moment, it appears that 
there is significant dispersion in dust radii as a function of age, with no strong correlation 
between these two quantities. 

5. Diverse Evolutionciry Histories? 

In summary, the dust masses of the 850 fim disks arc roughly consistent with the 
predictions of recent models for the coUisional evolution of planetesimal disks, in that the 
mass declines with time in the manner expected. However, the radii of the disks suggest 
a more diverse set of evolutionary histories for these systems that extend beyond histories 
dominated by self-stirring. 

Dominik & Decin (2003) similarly recognized the need for diversity in evolutionary 
histories in order to explain the large dispersion in the fractional luminosities {fd — Lm/ L^) 
observed by ISO in stars spanning a range of ages. In particular, in order to explain the 
large fractional luminosities [Lm/L^ — 10~^) observed by ISO in Vega-like systems both 
young (~ 10 Myr) and old (> 1 Gyr), they invoked the idea of "delayed stirring" in which 
a primoridal planetesimal disk is allowed to collide only after a specified waiting period. 
They speculated that since the Kenyon and Bromley models allow for the possibility of dust 
production over a long timescale, as protoplanet formation progresses from small to large 
radii, the models may describe the underlying physical mechanism that produces luminous 
dust disks spanning a wide range of ages. The distribution of dust disk radii as a function 
of age for the 850//m disks sample suggests that this is not the whole story. 

More specifically, some of the properties of individual systems are difficult to explain 
in detail solely with the self-stirring model. For example, HD104860 has the most massive 
(O.I6M0) as well as the largest dust disk. Compared to the ~ 40 — 50 AU disk sizes deduced 
for HD107146 and HD8907, the ~ 33 K dust temperature of HD104860 corresponds to a 
much larger orbital distance of ~ 110 AU. While the large disk size of HD104860 may reflect 
the more rapid collisional evolution (and inner disk clearing) expected for more massive 
disks, the large size cannot be accounted for with a simple scaling of the Kenyon & Bromley 
results. In their model, it would take ~ 400 Myr for the wave of collisional evolution to 
reach 100 AU even for a system with a mass 10 times the mass of the minimum mass solar 
nebula, an order of magnitude larger than the published age of 40 Myr for the system. 

These results suggest that processes beyond self-stirring play a significant role in the 
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evolutionary histories of planetesimal disks. One source of heterogeneity may be the influence 
of giant planets on planetesimal disks. Giant planets may clear regions of planetesimal disks 
and pump up the eccentricities and inclinations of planetesimals in other regions of disks, 
as well as sculpt the resulting dust debris. These effects have been invoked to explain the 
azimuthal asymmetries observed in the dust distributions surrounding some of the sources in 
the 850 iim disks sample (e.g., Vega; Wilner et al. 2002). It would not be surprising if giant 
planets were present in a large fraction of debris disk systems given the current detection 
statistics of precision radial velocity searches for extra-solar planets and the hkelihood that 
more giant planets remain to be discovered at larger orbital radii (Marcy et al. 2004). 

Giant planet formation is of course believed to have had a significant impact on the 
evolutionary history of our own Kuiper Belt. The fact that the dust content of the Kuiper 
Belt in the solar system is hkely located far below the points in Figure 4 at ~ 10~^Mq 
(Moro-Martin & Malhotra 2003a; Backman, Dasgupta & Stencel 1995) at an age of ~ 5 
Gyr, is likely the result of a more complex evolutionary history for the Kuiper Belt than 
the collisional grinding of a self-stirred planetesimal disk. Indeed, the low density of Kuiper 
Belt objects in the solar system today and the resulting low current rate of dust production 
is believed to be the consequence of giant planet formation which significantly depleted the 
planetesimal disk in the 30-50 AU region in the early solar system (e.g., Hahn & Malhotra 
1999). 

Since processes such as giant planet formation may produce significant dispersion about 
the evolutionary trend predicted by self-stirred models, large samples of debris disk systems 
spanning a large range in age may be needed to recognize the expected trend of dust radius 
increasing with age. This may be an interesting topic to return to when larger samples are 
available. 

6. Conclusions 

To conclude, we return to the question raised at the beginning of section 4.3, on the 
applicabihty of self-stirred models to the debris disks detected in our sample. As discussed 

above, there is some question as to whether self-stirred models apply. The sizes and ages 
of the dust disks surrounding HD8907 and HD107146 appear to fit the picture described by 
Kenyon & Bromley (2004). If the model is applicable to these systems, it suggests that the 
debris in these systems arises from massive planetesimal disks with initial masses several 
times the minimum mass solar nebula. In contrast, the radius of the dust disk surrounding 
HD 104860 appears too large to explain by a simple scaling of the model. More generally, 
the properties of the 850/im disks sample as a whole suggest the role of additional processes. 



- 21 - 



such as giant planet formation, in contributing to the production and shaping of the debris 
that is observed in at least some of the systems. 

On the one hand, this is encouraging for the FEPS program, since one of the aims 
of the program is to deduce the underlying giant planetary architectures of young solar 
analogues from the properties of their associated debris disks. On the other hand, the 
possibility that multiple processes affect the appearance of debris disks suggests that an 
important step toward the FEPS objective is to find ways to distinguish the relative roles of 
potentially multiple processes (e.g., embedded protoplanets, inner giant planets) in producing 
and shaping the debris. Spatially resolved observations of the debris may be one such 
discriminant. 

Our study made use of the significant work by the FEPS team that went into compiling 
and characterizing the FEPS sources. We would like to thank Steve Strom and Michael Liu 
for careful readings of the manuscript. JPW acknowledges support from NSF grant AST- 
0324328. This research has made use of the SIMBAD database and the Two Micron All 
Sky Survey, which is a joint project of the University of Massachusetts and IPAC/Caltech, 
funded by NASA and NSF. 
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TABLE 1 

Summary of JCMT/SCUBA measurements 



Source 


SpT 


Age 


d 


A 




Mb 






(Myr) 


(pc) 


(/xm) 


(mJy) 




IRXS J072343 


K3 


130 


24 


850 


< 3.9 


< 0.014 


HD 17925 


Kl 


100 


10.4 


850 


< 7.2 


< 0.005 


V383 Lac 


KO 


60 


50 


850 


< 8.1 


< 0.12 


HD 104860 


F8 


40 


47.9 


850 


6.8 ± 1.2 


0.16 










450 


47± 11 




HD 8907 


F8 


180 


34.2 


850 


4.8 ± 1.2 


0.036 










450 


22± 11 




HD 35850 


F7 


12 


26.8 


850 


< 5.4 


< 0.024 


HD 984 


F5 


40 


46.2 


850 


< 4.5 


< 0.059 


SAO 150676 


F5 


60 


78 


850 


< 6.3 


< 0.23 


HD 88638^ 


G5 


10 


37.5 


850 


< 4.8 


< 0.041 


HD 217343 


G3 


40 


32.0 


850 


< 6.0 


< 0.037 


HD 107146 


G2 


100 


28.5 


850 


20 ±3.2 


0.10 










450 


130 ± 12 




HD 166435'= 


GO 


10 


25.2 


850 


< 6.9 


< 0.027 


HD 77407 


GO 


30 


30.1 


850 


< 5.1 


< 0.028 



^ Upper limits are 3a, errors on the detections are ±a 

b K = 1.7 cm^ g-i; T = 50 K for upper limits, 33 K for HD 104860, 

48 K for HD 8907, and 51 K for HD 107146 

^ May be significantly older; dropped from FEPS sample subsequent 
to selection for this study. 
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TABLE 2 
(3 = 1 SED fits 



Source 




T 

dust 


Mdust 


Ldust 






(K) 


(Me) 




HD 104860 


best 


33 


0.16 


46 




min 


19 


0.34 


8 




max 


42 


0.12 


80 


HD 8907 




48 


0.036 


54 
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T I I 




A (yLim) 

Figure 1: Spectral energy distribution of HD 104860. The three filled squares are from the 
2MASS catalog. The star-disk system was not detected by IRAS; the open squares represent 
the FSC upper limits color-corrected for a dust temperature of 6000 K at 12 and 25 //m and 

40 K at 60 and 100 /xm. The filled circles show the SCUBA measurements at 450 //m and 
850 /im with ±a error bars. The solid line shows the best fit to the stellar photosphere and 
disk emission {Tdust = 33 K, /9 = 1), the dashed line shows the extrapolated contribution 
from the photosphere in the far-infrared. The grey area shows possible fits that do not 
violate the IRAS upper limits and lie within the systematic and calibration uncertainties of 
the SCUBA measurements. The corresponding range in dust mass and temperature is given 
in Table 2. 
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A (yLim) 

Figure 2: Spectral energy distribution of HD 8907. Detections are shown as filled symbols 
and include 2MASS photometry in the near-infrared, the 12 //m measurement from the IRAS 
FSC, 60 and 87 //m measurements from ISO, and the 850 //m SCUBA measurement. The 
2(7 SCUBA detection at 450 //m is shown as an open circle. For all points the error bars 
are ±a (when larger than the symbol size). The solid line shows the best fit to the stellar 
photosphere and disk emission {T^ust = 48 K, /3 = 1). 
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Figure 3: Spectral energy distributions for the 10 sources in the sample that were not 
detected in the submillimeter. The cluster of 3 points at 1 — 2 /im are fluxes from the 
2MASS catalog which were used to normalize Kurucz models appropriate for each stellar 
type. IRAS and ISO measurements from 12—100 fim are shown as solid squares for detections 
and open symbols with an arrow for upper limits. A few stars were observed at 1.3 mm by 
Carpenter et al. (2005), and these non-detections are also shown as open squares. The 3a 
SCUBA 850 iim upper limits from this work are shown as open circles. 
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Figure 4: Disk masses versus estimated stellar age for the sample of 13 objects discussed 
in this paper and for other sources in the literature. Detected sources are shown in black 
and grey points signify Sa upper limits based on a 50 K disk temperature, the average 
temperature for disks detected around solar type stars. The different symbols represent 
different stellar spectral types, as shown in the insert. The dashed line shows the slope 
expected if dust mass declines inversely with age and the darker shaded area shows the 
relatively older population in contrast to a significantly younger population in the unshaded 
region. The 70-300 Myr border between the two, due to low number statistics and age 
uncertainties, is shown by the lighter shaded vertical stripe. 
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Figure 5: Characteristic dust emission radii, assuming blackbody grains, plotted against 
estimated stellar age for the three detections in this paper and other sub-millimeter detections 
in the literature. The different symbols represent different stellar spectral types, as described 
in Figure 4. The upper limits represent the submillimeter-only excess sources from Wyatt 
et al. (2003). The trident symbol indicates the orbit of Neptune in the solar system. 



